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The Wiener-Khinchin theorem shows how the power spectrum of a stationary random signal 7( t) 
is related to its correlation function {I(t)I(t + t)). We consider non-stationary processes with the 
widely observed aging correlation function (7(t)7(t + r)} ~ t 7 </>EN(V/t) and relate it to the sample 
spectrum. We formulate two aging Wiener-Khinchin theorems relating the power spectrum to the 
time and ensemble averaged correlation functions, discussing briefly the advantages of each. When 
the scaling function </>en(*) exhibits a non-analytical behavior in the vicinity of its small argument 
we obtain aging 1/f type of spectrum. We demonstrate our results with three examples: blinking 
quantum dots, single hie diffusion and Brownian motion in a logarithmic potential, showing that 
our approach is valid for a wide range of physical mechanisms. 

PACS numbers: 05.40.-a,05.45.Tp 


Understanding how the strength of a signal is dis¬ 
tributed in the frequency domain, is central both in 
practical engineering problems and in Physics. In many 
applications a random process I ft) recorded in a time 
interval (0,t m ) is analyzed with the sample spectrum 
St m M = l/o m J 0)exp( — iu)t)dt\ 2 /t m , which is investi¬ 
gated in the limit of a long measurement time t m . For 
stationary processes, the fundamental Wiener-Khinchin 
theorem [l| relates between the power spectrum density 
and the correlation function C(t) = (I(t)I(t + r)} 

/•OO 

lim ( S tm (pj)) = 2 / C(t) cos(wr)dr. (1) 
tm ~>°0 J Q 


However in recent years there is growing interest in the 
spectral properties of non-stationary processes, where the 
theorem is not valid [M3- In general, there seems no 
point to discuss and classify spectral properties of all 
possible non-stationary processes. Luckily, a wide class 
of Physical systems and models exhibit a special type 
of correlation functions ( I(t)I(t + r)) ~ t" 1 (f^fr/t) for 
an observable lit) and the subscript EN denotes an en¬ 
semble average. Such correlation functions, describing 
what is referred to as physical aging, appear in a vast ar¬ 
ray of systems and models rangin g fr om glassy dynamics 


syst 

2, 11 jl3|, blinking quantum dots 14(, laser cooled atoms 


151 motion of a tracer particle in a crowded environment 
13 u. elastic models of fluctuating interfaces [18], de¬ 


terministic noisy Kuramoto models [191, granular gases 
, and deterministic intermittency }2ll | , to name only a 
few examples. In some cases the scaling function exhibits 
a second scaling exponent, (I(t)I(t + r)) ~ t 1 c/en(t /t^), 
or even a logarithmic time dependence [111] , however here 
we will avoid this zoo of exponents, and attain classifica¬ 
tion of the spectrum for the case /3 = 1. 

A natural problem is to relate between the sample 
spectrum of such processes and the underlying correla¬ 
tion function Q. That such a relation actually exists 
is obvious from the basic definition of the sample spec¬ 
trum, see Eq. m below. However, here we find a few 
interesting insights. First, the correlation function in its 


scaling form ( I{t)I{t + T )} ~ t 1 ^en(t / 1 ) is valid in Phys¬ 
ical situations, in the limit of large t and r. We here first 
formulate a theorem for ideal processes, where the aging 
correlation function is valid for all r and t , and then in the 
second part of the Letter, explore by comparison to re¬ 
alistic models the domain of validity of the ideal models. 
As a rule of thumb the aging Wiener-Khinchin theorem 
presented here for ideal models works well in the limit 
of low frequency. Further the limit of small frequency 
and measurement time t m being large is not interchange¬ 
able and should be taken with care. Secondly, the spec¬ 
trum in these processes depends on time t m , as already 
observed in I,®- The non-stationarity also implies a 
third theme, namely that the ensemble average correla¬ 
tion function is non identical to the time averaged correla¬ 
tion function, in contrast with the usual Wiener-Khinchin 
scenario. Thus we formulate two theorems, relating be¬ 
tween time and ensemble average correlation functions 
and the sample spectrum. The choice of theorem to be 
used in practice depends on the application. 

In physics the power spectrum is not only a measure of 
the strength of frequency modes in a system. Nyquist’s 
fluctuation dissipation theorem, for systems close to ther¬ 
mal equilibrium and hence stationary, states that the ra¬ 
tio between the power spectrum and the imaginary part 
of the response function, y(w), is given by temperature, 
i.e. ksT = itu)S(uj)/2^s [x(w)] (lljj. Similarly, effective 
temperatures are routinely defined by relating measure¬ 
ments of power spectrum and response functions of non- 
stationary processes [Ut^. Our goal here is to pro¬ 
vide the connection between the sample spectrum and 
the correlation functions, without which the meaning of 
the effective temperature becomes some what ambiguous. 
More practically, an experimentalist who uses the sample 
spectrum to estimate the spectrum of a non-stationary 
process, might wish to extract from it the time and/or 
the ensemble averaged correlation functions, and for that 
our work is valuable. 

Aging Wiener-Khinchin theorem for time averaged cor¬ 
relation functions. For a general process, the autocorre- 
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lation function (I(t)I(t + r)) is a function of its two vari¬ 
ables, unlike stationary processes, where the correlation 
function depends only on the time difference r. Using 
the definition of the sample spectrum we have 








( 2 ) 

We identify in this equation the ensemble average corre¬ 
lation function, but a formalism based on a time aver¬ 
age will turn out more connected to the original Wiener- 
Khinchin theorem, as we proceed to show. A change of 
variable t = fa — fa and relabeling integration variables 
gives 


find 


(S tm {u)) = 2(f m ) 1+T f df (1 — f) IfiTA {f) cos (wt m f). 
Jo 

( 6 ) 

This formula relates between the time average correla¬ 
tion function and the average of the sample spectrum, 
for ideal processes in the sense that we have assumed 
that the scaling of the correlation function holds for all 
times. It shows that the frequency uj times t m is the 
scaling variable of the power spectrum. 

Aging Wiener-Khinchin formula for the ensemble av¬ 
eraged correlation function. We now relate the power 
spectrum with the ensemble averaged correlation func¬ 
tion which has a scaling form 


= — / dT(t m -T){C TA (tm,r))cos(ujT). 
tm J o 

( 3 ) 

Here the time averaged correlation function is defined as 


CtA (trm T ) = 


tm T J o 


dfal(fa)/(fa+ t). ( 4 ) 


We now insert in Eq. ([3]) an aging correlation function 


(CTA(t m ,T)) = (im)'VTA (r/t m ) . 


(5) 


defining a new integration variable 0 < f = r /t m < 1 we 


{I(t + T)I(t)) = U^enCt/I)- (7) 


The two correlation functions are related with Eq. 
which upon averaging gives 

<?ta(s) = x 1 y{x) j ^ppdz ( 8 ) 

Jy(x) * +7 

with y{x) = x/{l — x). Considering the case 7 = 0 we 
insert Eq. (0 in Eq. 0 and find 


f 0 EA (U-l + 

Jo \t-x/ (us) 

(9) 

with w = u}t m . For the more general case 7 7 ^ 0 we show 
in the supplementary material that 


dec 




! + 1.1 2 + —• 
' 2 ’ 2 ’ ^ 2 ’ 



( 10 ) 


where 1 F 2 is a hypergeometric function and 7 > — 2 . 

This relation between the ensemble average correlation 
function and the sample averaged spectrum is useful for 
theoretical investigations, when a microscopical theory 
provides the ensemble average. Alternatively one may 
use the relation Eq. m to compute the time average 
correlation function from the ensemble average (if the 
latter is known) and then use the time averaged formal¬ 
ism Eq. 0 which is based on a simple cosine transform. 
The transformation Eq. m depends on 7 , which in ex¬ 
perimental situation might be unknown (though it could 
be estimated from data), while Eq. 0 does not, still 
both formalisms are clearly identical and useful. 

Relation with 1/f noise. We now consider a class of 
aging correlation functions, with the additional charac¬ 


teristic behavior for small variable r/t 


< mnt+T )) 




Aen — Hen 



( 11 ) 


Here Aen, Hen > 0,0 < v < 1,7 > — 1 and 7 — v > — 1 . 
Physical examples will soon follow. We use Eqs. 44070 
and by comparison of coefficients of small argument ex¬ 
pansion, we show that the time averaged correlation 
function has a similar expansion, with (Cta^j t)) ~ 

t ' 7 [A T a - B TA (T/t) 1 ' H-] with A T a = A E a/( 1+7) and 

Hta = Hea/(1 — ^ + 7 ). We can then insert this ex¬ 
pansion in Eq. 0 . by integration by parts and using 
uit m = 27t n 1 where n is an integer 


(s tm M> 


2 T (1 + v) sin ( 7 ^) H E a 
(7 - v + 1 ) (f m )‘'- 7 w 1 +" 


( 12 ) 


We see that the non-analytical expansion of the correla¬ 
tion function, in small argument, leads to a 1 // type of 
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noise, with an amplitude that depends on measurement 
time. Such an aging effect in the power spectrum was 
recently measured for blinking quantum dots 22 ], so this 
shall be our first example. 

Blinking quantum dots and trap model. As measure¬ 
ments show, blinking quantum dots, nano-wires and or¬ 
ganic molecules exhibit episodes of fluorescence inter- 
mittency, switching randomly between on and off states 


27 


. The on and off waiting times are random with 
a common power law waiting time distribution ip(r) ~ 
a behavior valid under certain conditions, like 
low temperature and weak external laser field. For this 
simple renewal model, and when the average on and off 
times diverge, namely 0 < a < 1 , we have 7 = 0 and the 
correlation function, with intensity in the on state taken 
to be Iq and in the off state to be zero, is 14] 


0en(z) = Iq 


1 sin(7ro;) 

2 4tt 


B 


\ + x 


; 1 — a, a 


(13) 


where x = r /t and B(z; a, b) is the incomplete beta func¬ 
tion. Importantly, this type of correlation function de¬ 
scribes not only blinking dots, but also the trap model, 
a well known model of glassy dynamics 0 . The con¬ 
nection between the two systems are the power law wait¬ 
ing times in micro-states of the system, though for the 
trap model a = T/T g where T is temperature, while 
0.5 < a < 0.8 in quantum dots experiments. Eq. m is 
valid only in a scaling limit for large r and t when the 
microscopic details of the model, e.g. the shape of the 
waiting time distribution ip{r) for short on and off blink¬ 
ing events, are irrelevant. Thus the scaling solution is 
controlled only by the parameter a. The time averaged 
correlation function is obtained from Eqs. (181131) 


0 ta(z)= 

Iq sin(7ra) 


T 2 
1 0 


(14) 


B (1 — x, a, 1 — a) 1 
1 — x a 


1 — x 


which is clearly non-identical to the corresponding en¬ 
semble averaged one. We may now use either Eq. © for 
the time average or Eq. ([9]) for the ensemble average to 
obtain St m (u>). Since 7 = 0 we use Eq. © and find 




sine 


(I) 


2 Q 


3 [M( 1 — a, 2; ioj)\ 


(!5) 

where M(a,b;z) is the Kummer confluent hypergeomet¬ 
ric function and 3 [.] refers to its imaginary part. We 
note that the sinc 2 (u)/ 2 ) term is the spectrum contri¬ 
bution from a constant (/ 2 ). As shown in Fig. |T] the 
spectrum Eq. m perfectly matches finite time simu¬ 
lation of the process where we used ip(r) = ar^( 1+Q ) 
for r > 1, a = 0.5, Iq = 1 and an average over 10 3 
on-off blinking processes was made. This indicates that 



FIG. 1. (color online). Power spectral density for the blinking 
quantum dot model with a = 0.5 with measurement time 
tm — 10 5 . Theory Eq. (1151) (red line) perfectly matches finite 
time simulation (blue o) and asymptotically the 1// noise Eq. 
(fl6l) (green line). 


the scaling approach works well, even for reasonable fi¬ 
nite measurement time. The theory predicts nicely not 
only the generic 1 // behavior but also the fine oscilla¬ 
tions and the crossover to the low frequency limit. As 
Fig. [T] demonstrates when ujt m is large, we get the 1 // 
noise result, which according to Eq. m is 

(StM) ~ (is) 

21 (1 + a) 

for 0 < a < 1. In this model v = 1 — a as the small argu¬ 
ment expansion of Eq. m shows. The asymptotic Eq. 
m agrees with previous approaches [ 3 } , the latter miss¬ 
ing the low frequency part of the spectrum, and the fine 
structure of the spectrum presented in Fig. [0 since for 
those non trivial aspects of the theory one needs the ag¬ 
ing Wiener-Khinchin approach developed here. Finally, 
we have assumed that the start of the blinking process 
at t = 0 (corresponding to the switching on of the laser 
field) is the moment in time where we start recording the 
power spectrum. If one waits a time t w before start of 
the measurement, the power spectrum will depend on the 
waiting time t w since the process is non-stationary Q ■ 
We note that a model with cutoffs on the aging behav¬ 
iors was investigated in [30j], in this case the asymptotic 
behavior is normal, namely the Wiener-Khinchin theo¬ 
rem holds. Indeed in experiments on blinking quantum 
dots with a measurement time of 1800 seconds the aging 
of the spectrum is still clearly visible [ 22 | , the latter mea¬ 
surement time is long in the sense that blinking events 
are observed already on the p, Sec. time scale. Hence 
cutoffs, while possibly important in some applications, 
are not relevant at least in this experiment. 
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Single file diffusion. We consider a tagged Brownian 
particle in an infinite unidimensional system, interact¬ 
ing with other identical particles through hard core col¬ 
lisions 31- 33|. This well known model of a particle in a 
crowded pore, is defined through the free diffusion coef¬ 
ficient D describing the motion of particles between col¬ 
lision events and the averaged distance between particles 
a. Initially at time t = 0 the particles are uniformly dis¬ 
tributed in space and the tagged particle is on the origin. 
In this many body problem, the motion of the tracer is 
sub-diffusive ( x 2 (t )) ~ a\fWpn\ft since the other parti¬ 
cles are slowing down the tracer particles via collisions 
31]. Normal diffusion is found only at very short times 
t < a 2 /D when the tracer particle has not yet collided 
with the other surrounding Brownian particles. Our ob¬ 
servable (so far called I(t)) is the position of the tracer 
particle in space x(t). The correlation function in the 

long time scaling limit is 0,113 


(x(t)x(t + t)) = a\l —\ft | 
7r 


,1 + i + 1 - 


(17) 


Such a correlation function describes also the coordinate 
of the Rouse chain model, a simple though popular model 
of Polymer dynamics [34]. Then by insertion and inte¬ 
gration we find using Eqs. m 


(tm) 3/2 \Jjfa2 ( S tm(u))= -z^Vd}{2 + COS{OJ)) (18) 

+ 2cos(w))C(v / 2w/tt) 

+ ^^' S (V / 2 d 7 V)(-<i + sin(w)), 

where the Fresnel functions are defined as C(u ) = 
/“ cos(nt 2 /2)dt and S(u) = /“ sin( 7 rf 2 / 2 )di. Generat¬ 
ing 10 3 trajectories of single file motion, for a system 
with 2001 particles, with the algorithm in H3 we have 
found the sample spectrum of the process x(t). As Fig. 
[ 2 ] demonstrates theory and simulation for (S tm (w)} per¬ 
fectly match without fitting. From the correlation func¬ 
tion Eq. m we have 7 = v = 1 /2 and hence according 
to Eq. (fl2l) 


<sm> - y^V 3/2 (19) 

for ujt m 1. This equation is the solid line presented in 
Fig. [3] which is seen to match the exact theory already for 
not too large values of u)t m - As in the previous example 
aging Wiener-Khinchin framework is useful in the pre¬ 
dictions of the deviations from the asymptotic result; as 
Fig. [2] clearly demonstrates some non-trivial wiggliness 
perfectly matching simulations. 

As mentioned in the introduction we have assumed a 
scaling form of the correlation function Eqs. (1517b . which 
works in the limit of t,r —> 00 . Information on the corre¬ 
lation function for short times is needed to estimate the 



FIG. 2. (color online). The power spectrum of tagged par¬ 
ticle motion x(t) with measurement times tm = 10 3 (blue •) 
and tm = 10 2 (pink x). For this single file model theory 
Eq. (1181) (red line) and asymptotic 1 // approximation Eq. 
nisi) nicely match simulation results. In simulation we used 
D = 1/2 and a = 1. Finite time deviations for tm = 10 2 are 
observed at high frequency as discussed in text. 


very high frequency limit of the spectrum. Hence the de¬ 
viations at high frequencies in Fig. 2 are expected. As 
measurement time is increased the spectrum plotted as 
a function of iot m perfectly approaches the predictions of 
our theory (see also the following example and Fig. 3). 

Diffusion in a logarithmic potential. While our pre¬ 
vious examples are based on long tailed trapping times 
and many body interactions which lead to a long term 
memory in the dynamics we will now briefly discuss a 
third mechanism using over damped Langevin dynamics 
in a system which attains thermal equilibrium. Con¬ 
sider the position x(t), which is the observable I(t), 
of a particle with mass in in a logarithmic potential 
U(x) = U 0 ln(l + x 2 )/2 


da; 

d t 


1 dU 
mj dx 


■V (*)■ 


( 20 ) 


Here the noise is white with mean equal zero satisfy¬ 
ing the fluctuation dissipation theorem and 7 is a fric¬ 
tion constant. Under such conditions the equilibrium 
probability density function is given by Boltzmann’s law 
P eq ( x) = exp[—U(x)/kbT]/Z , where Z is the partition 
function and T is the temperature. A key observation is 
that the potential is asymptotically weak in such a way 
that P eq ~ x ~u 0 /k b T £ or l ar ge x and for normalization to 
be finite we assume Uo/kbT > 1. The system thus ex¬ 
hibits large fluctuations in its amplitude in the sense that 
in equilibrium (x 2 ) diverges in the regime Uo/kbT < 3. 
Of course for any finite measurement time the variance 
of x(t), starting on the origin, is increasing with time 
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FIG. 3. (color online). Power spectrum of Brownian motion 
in a logarithmic potential (symbols) matches the asymptotic 
theory predicting 1// type of noise Eq. (1211) (green line). We 
use a = 1.75 ,t m = 5 • 10 2 (red A) and t m = 2 • 10 3 (blue 
o) with Uo = 1, 7 = 1 and m = 1. In addition we plot the 
theoretical zero frequency prediction (St. m ( 0)) (black line), see 
SM. The ensemble average was taken over 5000 realizations. 

but finite. Let a = Uo/{2kbT) + (1/2) and we focus on 
the case 1 < a < 2. The correlation function in this 
case was investigated in [35j. We here study only the 
1 // part of the spectrum, demonstrating the versatility 
of the theory using equation (fl2ll . since unlike previous 
cases the correlation function is cumbersome. To find the 
1 // noise we need to know, from the ensemble average 
correlation function, 7, v and Ben, while Aen must be 
finite (similar steps for other models will be published 
elsewhere [H|). As detailed in the SM [37] 7 = v = 2 — a 
and Ben = V^{^D) 2 ~ a ci/ [ZY(a)T{ 1 + a)] and Aen = 
Ben r(1 + a)/ [^tt(2 — a)ci] with D = kbT/mly the dif¬ 
fusion constant according to the Einstein relation, hence 
by using Eq. m we obtain 

( CX 7 T \ 1 

— J^EN^. (21) 

The constant Ci is given in terms of an integral of a spe¬ 
cial function [37|. We have simulated the Langevin equa¬ 
tion (USD and obtained finite time estimates for the power 
spectrum, which match the prediction Eq. m as shown 
in Fig. [3] without fitting. Importantly, the model of dif¬ 
fusion in logarithmic field is applicable in many systems, 
including diffusion of cold atoms in optical lattices [38|- 
Summary and Discussion. We have presented gen¬ 
eral relations between the sample spectrum and the 
time/ensemble averaged correlation function Eqs. (13191) . 
respectively. Those relations work for Physical models in 
the limit t m —> 00 and 10 —> 0 while the product ujt m is 
finite. In experiment t m might be long, but it is always 
finite. Hence, the theorem will work in practice in the 
low frequency regime. Indeed, a close look at Fig. [2] [3] 


shows finite time deviation at large frequencies, the ag¬ 
ing spectrum is approached when t m is increased. The 
fact that the scaled correlation function is observed in a 
great variety of different systems, serves as evidence of 
the universality of our main results, i.e. Eqs. USED- 
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